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ABSTRACT 


Application  of  the  low  level  wino  turbulence 

SPECTRUM  TO  THE  BALLISTIC  MISSILE  PROBLEM  WITH 

the  aid  of  Taylor's  hypothesis  is  discussed. 

Results  of  cross-spectral  analysis  are  discussed 
with  reference  to  predicting  the  wind  at  the 
launcher  from  a  sensor  at  SOME  point  in  space  away 
from  the  launcher.  Examples  of  spectral  and  cross- 
spectral  estimates  are  presented. 
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I  N  T  RODUCT  ION 


The  problem  of  obtaining  the  best  possible  measurement  of  the  low  level 

WIND  FOR  BALLISTIC  MISSILE  APPLICATION  HAS  NOT  BEEN  RESOLVED  TO  THE  SATIS¬ 
FACTION  OF  EITHER  THE  METEOROLOGIST  OR  THE  BALL  1  ST  I C I  AN .  THIS  WIND  PROBLEM 
IS  SEPARATED  INTO  TWO  AREAS,  THE  PPESHOOT  LAUNCHER  SETTING  FOR  WIND  EFFECT 
AND  THE  POST  LAUNCH  ANALYSIS  OF  THE  MISSILE  TRAJECTORY  AND  IMPACT. 

The  low  level  wind  measurements  are  taken  with  tower  mounted  sensors 

and/or  BALLOONS.  THESE  SENSORS  ARE  LOCATED  FROM  50  TO  2,000  FEET  FROM  THE 
MISSILE  LAUNCHER. 

Present  methods  cf  estimating  the  ballistic  wind  utilize  short  perioo 

(,  t..  INSTANTANEOUS  TO  ABOUT  ONE-MINUTE  AVERAGES)  WIND  DATA  PRIOR  TO  AND 
A1  IAUNCII  TIME.  MORE  ACCURATE  PRE-LAUNCH  ESTIMATES  OF  LAUNCH  TIME  W ' 

and  missile  dispersion  patterns  can  be  made  by  considering  more  stable  esti¬ 
mates  OF  MEAN  WIND  (I.E.,  10-  TO  60-M.NUTE  AVERAGES),  ^0  ALSO  ANALYZING 

WIND  DATA  TO  OBTAIN  SPECTRAL  AND  CROSS-SPECTRAL  CHARACTERISTICS  OF 
TURBULENT  WIND  FIELD. 

DISCUSSION  .  _ 

- PuLER I  AN  TIME  CORRELATION  CONSIDERS  VELOCITY  FLUCTUATIONS  AT  A  FIXED 

POINT  WITH  RESPECT  TO  TIME.  EULERIAN  SPACE  CORRELATION  CONSIDERS  SIMULTA¬ 
NEOUS  VALUES  OF  VELOCITY  FLUCTUATIONS  AT  TWO  OR  MORE  POINTS  IN  SPACE  WITH 
RESPECT  TO  THE  DISTANCE  BETWEEN  THE  POINTS. 

Nearly  all  wind  measurements  are  of  the  Euler ian-t ime  form,  but  wind 

FLUCTUATIONS,  AS  THEY  AFFECT  MISSILES  IN  FLIGHT,  ARE  OF  THE  EULER  I AN-SPACE 
FORM,  I TURBULENCE  DESCRIBED  IN  FORMS  OF  WAVE  LENGTH  OR  WAVE  NUMBER  IN 
AN  EFFECTIVELY  FROZEN  ATMOSPHERE.  To  RESOLVE  THIS  WE  USE  G.  I.  TAYLOR  S 
HYPOTHESIS  (Ref.  1),  WHICH  TRANSFORMS  A  TIME  SERIES  INTO  A  SPACE  SERIES. 

The  hypothesis  is:  "If  u,  the  component  at  a  fixed  point  of  turbulence 
motion  in  the  direction  of  the  main  stream  in  a  wind  tunnel,  is  resolved 
into  harmonic  components, the  mean  value  of  u2  may  be  regarded  as  being  the 
sum  OF  contributions  from  all  frequencies.  If  u2F(n)  on  is  the  contribu¬ 
tion  from  frequencies  between  n  and  n  +  cn,  then 


F(n)  on  =  1 


(1) 


"If  F(n)  is  plotted  against  n,  the  diagram  so  produced  is  a  form  of 
the  spectrum  curve. - 


1 


"It  is  clear  that  when  the  eddies  are  large  the  correlation  Rx  between 

SIMULTANEOUS  VALUES  OF  U  AT  DISTANCE  X  APART  MUST  FALL  AWAY  WITH  INCREASING 
X  MORE  SLOWLY  THAN  WHEN  EDDIES  ARE  SMALL.  ONE  MAY  THEREFORE  ANTICIPATE 
THAT  WHEN  THE  (Rx,  x)  CURVE  HAS  A  SMALL  SPREAD  IN  THE  X  CO-ORDINATE  THE  F(*l) 
CURVE  WILL  EXTEND  TO  LARGE  VALUES  OF  n  OR  VICE  VERSA. 

"If  THE  VELOCITY  OF  THE  AIR  STREAM  WHICH  CARRIES  THE  EDDIES  IS  VERY 
MUCH  GREATER  THAN  THE  TURBULENT  VELOCITY,  ONE  MAY  ASSUME  THAT  THE  SEQUENCE 
OF  CHANGES  IN  U  AT  THE  FIXED  POINT  ARE  SIMPLY  DUE  TO  THE  PASSAGE  OF  AN 
UNCHANGING  PATTERN  OF  TURBULENT  MOTION  OVER  THE  POINT,  I.E.,  ONE  MAY  ASSUME 
THAT 


U  =  ♦(T)  -  ,  (7) 

WHERE  X  IS  MEASURED  UPSTREAM  AT  TIME  T  =  0  FROM  THE  FIXED  POINT  WHERE  U  IS 

MEASURED.  I'N  THE  LIMIT  WHEN  u/l)  ^.0  (j)  IS  CERTAINLY  TRUE.  ASSUMING 

THAT  (7)  IS  STILL  TRUE  WHEN  u/U  IS  SMALL  BUT  NOT  ZERO,  Rx  IS  DEFINED  AS 


Rx 


(8)" 


U  =  MEAN  WIND  STREAM  VELOCITY. 

From  this  is  derived  the  relationship 


/  _x_ 

Rx  =  »  (J)  »  VT  +  U 


F(n)  cos  — do 


AND  USING  THE  FOURIER  INTEGRAL  THEOREM  THE  FOLLOWING  EQUATION  IS  FOUND  FOR 

F(n) 


F(n ) 


)  r  00 

=  -  /  Rx  cos  g  n  n  X  DX  . 

u  Jo  U 


Thus  if  F(n)  is  known  Rx  may  be  calculated  or  vice  versa. 

Taylor's  hypothesis  has  been  proven  experimentally  in  wind  tunnels. 


2 


Other  experiments  have  demonstrated  that  the  hypothesis  x  =  Ut  holds  true 
IN  the  atmosphere  for  wavelengths  of  less  than  1,000  FEET  ( Re F .  2,  3,  4,  5)' 
Sufficient  experiments  have  not  been  made  to  verify  the  hypothesis  for 

WAVELENGTHS  IN  EXCESS  OF  1,000  FEET  IN  THE  ATMOSPHERE. 


POWER  SPECTRAL  ESTIMATES 

0  F 

TURBULENCE  STRUCTURE 


In  recent  years  considerable  work  has  been  done  on  THE  evaluation  of 
the  power  spectral  estimates  of  the  turbulent  wind  field  at  a  single  point. 
Since  this  can  be  transformed  into  space  spectral  estimates  by  Taylor’s 
hypothesis^  the  spectral  estimates  of  time  series  shall  be  DISCUSSED  r I rst . 
In  all  following  examples  the  CRITERION  of  u/u  being  small  is  approximately 
MET,  I  .E  .,  GUST  I  NESS  RATIO  (  CT  U/U)  <  O.3. 

Huss  and  Bushnell  (Ref.  6)  calculated  the  power  spectrum  " -  by 

FINDING  THE  VARIANCE  OF  VELOCITY  AS  A  FUNCTION  OF  I  AND  J,  TWO 

LENGTHS  OF  TIME  INTERVALS,  AND  THEN  DIVIDING  BY  THE  CORRESPONDING  BANDWIDTH 

1  -  1  I 

'g  ,  “  . ~2j  I  AND  FINALLY  ASSIGNING  THE  RESULT  TO  THE  MIDFREQUENCY 


1  r  j _ +  j _ 

2  L  2 i  2j  _ 


WERE  FOUND  FROM  THE  RELATION 

2  ? 
v  t  /  \  -  cr  ,  ,  .  • 

('/'-)  (  j/l  ) 

L  =  FUNDAMENTAL  SAMPLE  PERIOD. 

Spectral  estimates  may  be  normalized  by  dividing  by  the  square  of  the 

MEAN  WIND  SPEEO. 

Figures  1  and  2  are  spectral  estimates  of  the  wind  15  feet  and  180  feet 
above  the  surface  * t  White  Sands  Missile  Range  computed  using  the  method 
described  by  Huss  „..j  Bushnell. 


The  variances 
CT2  («/o)  = 


3 


5  6  7  B  9  1*-vrw  2  3  4  567891 

Period  of  Eddy  (Sec.) 


3  4  5  6 


Fig.  1 
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Power 

COVAR I ANCE 


spectral  estimates  F(n)  are 

FUNCTION  (Rk): 


ALSO  OBTAINED  BY  USE  OF  THE  AUTO- 


R 


K 


M  -  K 

y  xx-  '  j 

L  *  |+K 

i  =  i 


and  the  Fourier  transformation  (Ref.  8,  9): 


F 


(n 


o) 


(R  +  Rm>  + 

o  M 


M  -  1 

R„  +  R»]*  Z  R«  C0S 

K  =  1 


M  -  1 


F(rV~ 


(-1)M  Rk 


i 

J 


£  (-i)K  rk'| 


WMCRE 


K  IS  THE  LAG  NUMBER, 

M  IS  THE  MAXIMUM  LAG  NUMBER  USED, 
N  IS  THE  SAMPLE  SIZE, 


6 


X  IS  THE  INDIVIDUAL  DEVIATION  FROM  A  MEAN, 
n  IS  FREQUENCY  DETERMINED  BY  CHOICE  OF  K. 
The  bandwidth  is  given  by  the  equation 


An  .  -(yjjgJ 

2  AtM 

AND  THE  CENTER  FREQUENCY  IS  GIVEN  BY 


WHERE  A  T  IS  THE  OBSERVATIONAL  INTERVAL.  THE  SPECTRAL  ESTIMATES  ARE  THEN 
SMOOTHED  BY  A  WEIGHTED  MOVING  AVERAGE. 

Panofsky  and  Deland  (Ref.  10)  found  from  spectral  analysis  of  wind 
DATA  AT  BROOKHAVEN,  LONG  ISLAND,  AND  O'NEILL,  NEBRASKA,  THAT,  "THE  SPECTRUM 
OF  LATERAL  VELOCITY  COMPONENTS  CAN  BE  DIVIDED  MOST  CLEARLY  INTO  LOW-FREQUENCY 
CONVECTIVE  AND  HIGH-FREQUENCY  MECHANICAL  PORTIONS.  THE  CONVECTIVE  PORTION 
IS  ALMOST  ENTIRELY  A  FUNCTION  OF  LAPSE  RATE  OR  SHORT-WAVE  RADIATION,  WITH 
A  TENDENCY  TO  INCREASE  WITH  HEIGHT.  It  IS  ESSENTIALLY  INDEPENDENT  OF  WIND 
SPEED  AND  GROUND  ROUGHNESS.  THE  MECHANICAL  PORTION,  ON  THE  OTHER  HAND,  IS 
SENSITIVE  TO  GROUND  ROUGHNESS  AND  INDEPENDENT  OF  STABILITY,  AND  TENDS  TO 
DECREASE  WITH  HEIGHT.  SINCE  THE  CONVECTIVE  PART  OF  THE  SPECTRUM  CAN  BE 
LARGE  AT  DAYTIME,  THE  TOTAL  VARIANCE  OF  LATERAL  VELOCITY  SHOWS  A  TREMENDOUS 
DIURNAL  VARIATION. 

"The  PROPERTIES  OF  THE  SPECTRUM  OF  THE  LONGITUDINAL  WIND  COMPONENT 
ARE  SIMILAR  TO  THOSE  OF  THE  LATERAL  COMPONENT.  HOWEVER,  THE  LOW-FREQUENCY 
PORTION  OF  THE  SPECTRUM  IS  CONSIDERABLE  EVEN  iN  STABLE  AIR,  SHOWING  THAT 
THE  LARGEST  EDDIES  AT  NIGHT  ARE  ELONGATED  ALONG  THE  WIND."  FIGURE  3  IS  A 
SCHEMATIC  PRESENTATION  OF  SPECTRA  AT  O'NEILL.  CRAMER  (REF.  11)  HAS  FOUND 
THAT  THE  MECHANICAL  TURBULENCE,  I  . E .,  FRE QUENC  I ES  HIGHER  THAN  1  Cy/mIN, 

DEPENDS  PRIMARILY  ON  THE  SQUARE  OF  THE  MEAN  WIND  SPEED.  CONVECTIVE  TURBUL- 
LENCE,  I  ,E .,  FREQUENC IES  LOWER  THAN  1  Cy/m'N,  DEPENDS  UPON  THE  MEAN  WIND 
SPEED  AND  THERMAL  STRATIFICATION.  He  ALSO  FOUND,  AS  THEY  DID  AT  BROOKHAVEN, 

(Ref.  12)  that  the  magnitude  or  standard  deviation  of  the  fluctuations  of 

THE  WIND  DIRECTION  IS  A  GOOD  INDICATOR  OF  THERMAL  STRATIFICATION  AND,  THERE¬ 
FORE,  IS  A  GOOD  PREDICTOR  OF  CONVECTIVE  TURBULENCE  SPECTRA. 

IN  THE  ATMOSPHERIC  TURBULENCE  THERE  IS  NO  LOWER  FREQUENCY  LIMIT  SO 
PRACTICAL  LOW  FREQUENCY  LIMITS  MUST  BE  SELECTED.  A  LOGICAL  CUTOFF  VALUE 
WOULD  BE  BETWEEN  THE  CONVECTIVE  TURBULENCE  AND  THE  DIURNAL  VARIATIONS  WHERE 
THERE  IS  A  MARKED  DECREASE  IN  THE  SPECTRAL  ENERGY.  THIS  DECREASE,  OR  GAP, 


8 


OCCURS  IN  THE  FREQUENCY  RANGE  OF 

Figure  4  (Ref.  13)»  Frequencies 

ENERGY  WILL  NOT  BE  CONSIDERED  IN 


1  TO  2  CYCLES  PER  HOUR  AS  INDICATED  IN 
LOWER  THAN  THOSE  PRODUCED  FROM  CONVECTIVE 
THIS  PAPER. 


CORRELATION  OF  WIND 
BETWEEN  TWO  POINTS 

- —  .  —I.  ....  . . . .  .in 


Cross-spectral  analysis  determines  the  correlation  of  various  frequen¬ 
cies  BETWEEN  TWO  TIME  SERIES  AND  IS  COMPOSED  OF  TWO  PARTS,  THE  COSPECTRUM 
AND  THE  QUADRATURE  SPECTRUM.  THE  COSPECTRUM  PORTION  CONSIDERS  THE  SIMULTA¬ 
NEOUS  RELATIONS  BETWEEN  TWO  SERIES  AND  THE  QUADRATURE  SPECTRUM  CONSIDERS 
THE  RELATIONSHIP  BETWEEN  TWO  SERIES  WHEN  ONE  IS  OUT  OF  PHASE  BY  A  QUARTER 
PER IOD. 

The  COHERENCE  FUNCTION,  WHICH  IS  ANALOGOUS  TO  THE  SQUARE  OF  THE  COR¬ 
RELATION  COEFFICIENT,  EVALUATES  THE  RELATIONSHIP  BETWEEN  TWO  TIME  SERIES 
FOR  VARIOUS  FREQUENCIES  AND  IS  DEFINED  BY 


COH 


(n) 


(n) 


(n) 


X/Y  F(n)vF(n)< 


WHERE 


C(n)  ls  THE  COSPECTRUM  BETWEEN  THE  TWO  TIME  SERIES  X  AND  Y  AND 


Q,  \ 

(n)  IS  THE  QUADRATURE 
COMPUTATIONAL  PROCEDURES  OF 
10. 


spectrum  between  the  same  two  series.  For 

THE  CROSS  SPECTRUM  SEE  REFERENCES  5>  8,  9 >  °R 


BROOKHAVEN  HAS  APPLIED  THE  C R OSS-S PE CTR AL  EQUATIONS  TO  WIND  PROFILE 
PREDICTION  TECHNIQUES  ( R,E  F .  1 4—  1 7 )  USING  THE  75"  >  150-  AND  300-FOOT  LEVELS 

ON  ONE  TOWER  (Ace)  AND  THE  150-FOOT  LEVEL  ON  A  SECOND  TOWER  ( K I NG )  LOCATED 
900  FEET  AWAY  FROM  THE  FIRST.  TABLES  I  AND  II  ARE  COHERENCE  VALUES  COMPUTED 
AT  Brookhavfn  (Ref.  15)  where  the  wind  was  blowing  directly  from  one  tower 

TO  THE  OTHER.  BOTH  RUNS  WERE  UNDER  LAPSE  CONDITIONS  WITH  THE  AZIMUTH 
FLUCTUATION  RANGING  FROM  15°  TO  45°. 

(Cont.  on  Page  13) 


9 


B— M— i—MBB—l 


3eBBaaa«»ii _ 

—  — —  MMIf  j 

wBayjMil 

gsaasrissaf  I 


SiSSsHilaEMsS^g^agsggg^fiagiiiigggSsS^yfc-sggsiiiisgiiriMml 


MwoniKBrflBn 

Safiai-  jcror 


■■  swse  mafias  bbosi  53f?---!3r-aKfi  he 

ssi  Ssfeisfe  £12*31  ar.-SSII-CaiS  | 

^fBB—giBBB55g  aggespaagB  gsgugsagsa  asggg  Eggsa  Fcsar  ~  ■ 

iMHtM  sesiMaec  iM  isi2§isii§ — »  ^^gggisi 


IIHilBHIL:lBI^BB^fcs:sss==3Bassssss====s=ss=srr-sSssss5ssss=s=sik: 


*r  '-’Ut  c 

IMMP  WBHB  —UBS 

■  IjESHr.mar’  ~=r 


:sccdsiR 


I  P-OBb  aUkae  taem  giiEi  lanes  esWl 

iMMMmnr  aui«HaHnr 
HE’iSHS- SBEBr  - aBf -Hi  L’E ■ 


lauaLE  JC  JC.C- 

■mbse  bhb  sb  amt  nra  a  asms  amssas  SSBB  m  bpcbc  c 
ii-j— —— — — ■  -— .^b—w—  BgiaaiB— gacqcpgg^g 
ii<imr  TasawiB^Mw— ■■— iaw BiHBgg^sgj 

i§iiiii!!pg5iiiipp  pin  5=111  ^n== g==== i§=si?i=p  mu  ===51 1 

j  5|==||ii=| =f===Ei=55  s==i§  =====  === 


siSSSSs 


ji  ms usi 


gesjjg  agsijj 


I— — BW 


JMHHIIHI 


is  ggn^g  aggBBaaHBa 


|ME= 


■MMHianii 


M 


2=55=5 


nmsanBi 


sail 


:========: 


nmj 


3  Secajod  assess 

MHHHnnii 


mEHiMm  min 


aasra 


jtiw  ats£i=y  asaaas 


:gnm  esssa  sen  nr  gnRn&r-v.nw  asn&s 


i  on 


:sssl 


HIS 


net 


iaEIBMI 


}£9ifE1k!B££ 


■■■I  ■■■—■IIM II  ■■■■■■!■!  I 


Ifl  BBE3H  Ej 

II  »£3l  ram  men  grass  r 


in 


las 


::===: 


— — Ml 

——I 

iil—IMIIinga 


1  IT  ■■■■-■■■  MW 


HHHE 


5SB5B3- 


■MIMilMMMCI 


nrai 


1  mMmmm  in§i  urn  ii=si  mu  i 


■■■■■■ ■■■■■  ru 


IMM1UI 


■■mi  — im  aggafettBygg  sessa  bsissb  granges 

limM 

|ii£gEl 

HHys9 

IliHIiiiiillllii 


bhi 


-saasatl 


:  =  =  ■  :  . - 


6  g^aS/gM 6 


TABLE  I 


COHERENCE  VALUES  - 

HIGHER  WIND  SPEEDS 

U  150*  FT 

=  6.55  m/sec, U  150 

FT  =  6.02  M/Sec 

:,U  300  ft  =  6.58  M/S 

Coherence 

(After 

Singer) 

Freq  (cy/hr) 

150*  -  150 

150*  -  300 

150  -  300  (ft) 

0 

0.25937 

O.1536I 

0.18454 

5 

0.77970 

o.64o64 

0.75488 

10 

0.46451 

0.57854 

0.56491 

15 

0.38664 

0.32281 

0.38890 

20 

0.19913 

0.18186 

0.24469 

25 

0.199^ 

0.20020 

0.32738 

30 

0.51674 

0.21136 

0.33861 

35 

o.4 1642 

0.32838 

0.38582 

4o 

0.25262 

o.4o498 

0.27840 

45 

0.09653 

0.19290 

0.20341 

50 

0.06619 

0.08499 

0.20543 

55 

0.07004 

0.07251 

0.12810 

60 

0.15204 

0.07101 

0.18397 

80 

0.13694 

0.03565 

o.o4i97 

100 

o.o4i59 

0.05613 

0.00115 

120 

0.16702 

0.05615 

0.03449 

l4o 

0.07087 

0.04435 

0.07055 

160 

0.04970 

0.02469 

0.00059 

180 

0.02866 

0.07614 

0.08271 

*  King  Tower. 
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TABU  I  I 


COHERENCE  VALUES  -  LOWER  WIND  SPEEDS 
U  150*  FT  =  3.71  M/SEC^U  150  FT  =  3.67  M/SEC^U  300  FT  =  3. 70  M/SEC 


Coherence  (After  Singer) 

Freq  (cy/hr)  150*  -  150  150*  -  300  150  -  300  (ft) 


0 

0.32445 

0.18109 

0.27907 

5 

0.73456 

0.6l66l 

0.78931 

10 

O.57689 

0.44867 

0.51747 

15 

0.20045 

0.04538 

0.26880 

20 

0.22248 

0.06281 

0.27279 

25 

0.l6l6l 

0.07387 

o.o4oo9 

30 

O.06685 

o.o4o97 

0.00957 

35 

0.00713 

0.15390 

0.00204 

4o 

0.04656 

0.12309 

o.oo4o4 

t 

45 

0.17881 

0.11434 

0.00005 

50 

O.I6562 

0.07509 

0.04376 

55  ' 

0.03171 

0.07824 

0.02927 

60 

0.05403 

0.08705 

0.01995 

80 

0.03207 

0.11220 

0.12661 

100 

0.05224 

0.14466 

0.09286 

120 

0.11845 

0.09259 

0.17874 

i4o 

O.09658 

0.06787 

0.12745 

160 

0.01974 

0.02632 

0.11042 

180 

0.00346 

0. 

0.00905 

*  King  Tower 
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TABLE  I  I  I,  ALSO  EXTRACTED  FROM  REFERENCE  1  5>  CONSISTS  OF  COHERENCE  VALUES 
BETWEEN  THREE  LEVELS  ON  A  SINGLE  TOWER  UNDER  A  NIGHTTIME  (TEMPERATURE 
INVERSION)  CONDITION.  THE  DIRECTION  TRACE  FOR  THIS  RUN  APPROXIMATES  A 
SMOOTH  CURVE. 

Cramer  (Ref.  5)  analyzed  the  data  obtained  from  Project  Prairie  Grass 

FOR  CORRELATIONS  OF  THE  U  (DOWNWIND)  AND  V  (CROSSWIND)  FLUCTUATIONS  OF  THE 
WINO  AT  A  HEIGHT  OF  TWO  METERS.  He  DETERMINED  THAT  THE  QUADRATURE  SPECTRAL 
ESTIMATES  WERE  NOT  STATISTICALLY  SIGNIFICANT  AND  SO  COMPUTED  THE  CORRELA¬ 
TIONS  (Rq)  BY  THE  EQUATION 


RC(n)  " 


Figure  5  (Cramer,  Ref.  5)  presents  correlation  coefficients  vs  separation 
DISTANCE  FOR  THE  VARIOUS  PERIODS  GIVEN  IN  TABLE  IV.  THIS  FIGURE  REPRESENTS 
SAMPLES  WHERE  THE  WIND  WAS  BLOWING  PARALLEL  TO  THE  SENSOR  ARRAY. 

Coherence  values  between  various  levels  on  the  White  Sands  meteorologi¬ 
cal  TOWER  DURING  DAYTIME  (LAPSE)  CONDITIONS  ARE  PRESENTED  IN  TABLE  V.  THE 
POWER  SPECTRA  FOR  THIS  RUN  ARE  GIVEN  IN  FIGURES  1  AND  2. 

Some  other  approaches  to  relate  wind  measurements  at  one  point  to  those 

AT  ANOTHER  POINT  IN  SPACE  ARE  DISCUSSED  BY:  HOWCRAFT  AND  SMITH  (f?EF.  l8); 

Cramer,  Record,  and  Vaughn  (Ref.  19);  Lamberth  and  VeitH  (Ref.  20);  Durst 
(Ref.  21);  and  Court  (Ref.  22). 


APPLICATION  OF  WIND  SPECTRA 
INFORMATION  TO  BALLISTICS 


Zbrozek  (Ref.  23)  and  Thorson  and  Bohne  (Ref.  24)  discuss  the  output 

POWER  SPECTRUM  OF  AN  AIRCRAFT  OR  MISSILE  [^(n)]  AS  A  FUNCTION  OF  THE 
POWER  SPECTRUM  OF  TURBULENCE  [F(fi)  ]  AND  THE  AIRCRAFT  FREQUENCY  RESPONSE 
FUNCT  ION  [  T(R)  ]  . 

The  OUTPUT  SPECTRUM  CAN  BE  CALCULATED  BY  the  RELATIONSHIP 

Fm  («)  =  [T(«)]2  x  F(fi) 
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(Cont  on  Page  17) 


TABLE  I  I  I 


COHERENCE  VALUES  -  TEMPERATURE  INVERSION  CONDITIONS 
U  75  FT  =  3.44  M/Sec,  U  150  FT  =  6.33  M/SEC, U  300  FT  =  9. 01  M/SEC 
Coherence  (After  Singer) 


Freq  (cy/hr) 

75  -  150 

75.  -J300 

150  -  300  (ft) 

0 

0.44598 

0.13006 

0.14846 

5 

0.60920 

0.17281 

0.31483 

10 

0.37815 

0.35140 

0.35924 

15 

0.17952 

0.09043 

0.15680 

20 

0.28390 

0.06399 

0.14721 

25 

0.13622 

0.02686 

0.07647 

30 

0.07008 

0.00058 

0.12540 

35 

o.i4oo3 

0.00111 

0.12700 

4o 

0.13978 

o.oo489 

0.11364 

45 

0.09225 

0.02359 

0.05435 

50 

0.24533 

0.12522 

0.00106 

55 

0.08598 

0.04293 

0.15396 

60 

0.03975 

0.0451 4 

o.oi4o4 

80 

0.00212 

0.02632 

0.07168 

100 

0.10833 

0.11111 

0.09259 

120 

0.06974 

0.06579 

0.05625 

i4o 

0.03137 

0.02941 

0.10667 

160 

0.02564 

o.o64io 

0.02778 

180 

0.01070 

0.02941 

0.09191 

TABLE  IV 


CENTRAL  FREQUENCIES  ECR  SELECTED  LAG  VALUES 
k  Tr  =  1 —  (sec) 


1  128 

2  64 


3 

4 

5 

6 
8 


43 

32 

26 

21 

16 


10 


12.8 


12 


10.7 


15  8.5 

20  6.4 


24 


5-3 


Central  frequencies  of  frequency  intervals  associated  with  selected 

VALUES  OF  K  ( LAG  )  USED  IN  OBTAINING  SCALE  ESTIMATES;  FOR  CONVENIENCE, 
DATA  ARE  INVERTED  AND  EXPRESSED  IN  TERMS  OF  PERIOD  RATHER  THAN  FREQUENCY. 
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COHERENCE 


k  Period  (sec) 


1 

24o 

2 

120 

3 

8o 

4 

6o 

5 

48 

6 

4o 

8 

30 

10 

24 

12 

20 

16 

15 

1 

240 

2 

120 

3 

8o 

4 

6o 

5 

48 

6 

4o 

8 

30 

10 

24 

12 

20 

16 

15 

TABLE  V 


VALUES 

DURING  LAPSE 

CONDITIONS  - 

WSMR 

Coherence  (u  component) 

15  -  63  15  -  111  15  -  203 

63  -  HI 

63  -  203  (ft) 

.7006 

•733*1 

.25*19 

.9173 

.7066 

.5871 

.3781 

.1146 

.7517 

•5531 

.7257 

.3762 

.3522 

.6187 

•3540 

.7^3 

.2596 

.5112 

.4579 

.5177 

.3307 

.0066 

.1585 

•1559 

.3720 

.2991 

.0355 

.0300 

.0914 

•0553 

•0533 

.5162 

.4415 

.3262 

.0066 

.0308 

.5389 

.0423 

.0029 

.3072 

.2638 

.2560 

.1863 

.0962 

.2250 

.0053  .2938  .1622 

Coherence  (v  component) 

.0260 

.1*137 

.2200 

.15*17 

.2064 

.9809 

.9945 

.0945 

.0999 

.0271 

.8032 

-=*- 

ON 

1 — 

• 

.0281 

,o468 

.1488 

.3229 

.0318 

.2019 

,i4io 

.1107 

.5752 

.0853 

.0801 

.0230 

.4179 

.3273 

.0128 

.0586 

.0900 

.6101 

.2747 

.0097 

.0001 

.2164 

•3395 

.1626 

.1340 

.1214 

.1675 

.2121 

.3859 

•  1157 

.1556 

.1995 

.0229 

.0335 

.2778 

.0438 

.0562 

.1303 

.0655 

.4364 
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WHERE 


ft  =  SPACE  FREQUENCY  OR  WAVE  NUMBER.  FOR  WIND  FIELD  CHARACTERISTICS 
ALONG  A  HORIZONTAL  PATH  [F(fi),fi]  CURVES  MAY  BE  ESTIMATED  USING  TAYLOR'S 
HYPOTHESIS  AND  AVAILABLE  f(n)>  11  ]  CURVES. 

Approximate  [F(ft)>  3  curves  for  a  vertical  path  may  be  directly 

ESTIMATED  FROM  PILOT  BALLOON  WIND  MEASUREMENTS. 

When  the  admittance  function,  i .e .,  response  to  a  particular  space 
FREQUENCY,  OF  A  ROCKET  IS  KNOWN  WE  CAN  DETERMINE  THE  BANDWIDTH  OF  THE 
TURBULENCE  SPECTRUM  THAT  WILL  CAUSE  THE  MISSILE  TO  DEVIATE  FROM  ITS 
EXPECTEO  TRAJECTORY.  FOR  EXAMPLE,  SUPPOSE  AN  OVERDAMPED  MISSILE  WITH  A 
NATURAL  WAVELENGTH  OF  1,000  FEET  IS  FIRED  HORIZONTALLY  AND  HAS  AN  ADMIT¬ 
TANCE  FUNCTION  DESCRIBED  BY  FIGURE  6.  USING  THE  TURBULENCE  SPECTRAL  VALUES 

in  Figure  1  with  the  assumption  that  Taylor’s  hypothesis  holds  true,  and 

THE  ADMITTANCE  RATIO  VALUES  IN  FIGURE  6  WE  OBTAIN  AN  EFFECTIVE  SPECTRAL 
CURVE  [  T(ft),  ft  ]  AS  SEEN  IN  FIGURE  7*  Fr0M  TH,S  CURVE  WE  CAN  DETERMINE 
THE  HIGHEST  FREQUENCfES  THAT  AFFECT  THE  MISSILE  DISPERSION  PATTERN. 

Smooth 'ng  and  filtering  techniques  may  be  employed  to  eliminate  the  fre¬ 
quencies  WHERE  THE  ENERGY  LEVEL  IS  SO  LOW  THAT  IT  IS  INSIGNIFICANT. 

Smoothing  and  filtering  of  data  is  discussed  by  Holloway  (Ref.  2F). 

Space  spectral  information,  when  introduced  into  ballistic  equations, 

WILL  DESCRIBE  THE  EFFECT  OF  LOW  LEVEL  WIND  ON  THE  DISPERSION  PATTERN  OF 
MISSILE  TRAJECTORIES.  BALLISTIC  EQUATIONS,  USING  THE  TURBULENCE  SPECTRA, 

ARE  BEING  DEVELOPED  BY  RaCHELE  ( Re F .  26)  AND  WALTER  (REF.  27). 

When  a  missile  is  to  be  lauw-  ,  low  level  wind  is  measured  at  a 

POINT  OR  POINTS  IN  SPACE  SOME  D  .  E  FROM  THE  LAUNCHER  AND  THE  ANTICI¬ 

PATED  MISSILE  PATH.  THIS  WIND  !  M*-  ,'RMA  T  )  CN  IS  USED  TO  CORRECT  THE  LAUNCHER 
AZIMUTH  AND  ELEVATION  FOR  THE  WIND  EFFECTS  AND  IS  ALSO  USED  FOR  POST  FIRING 
ANALYSIS  OF  THE  FLIGHT. 

A  METHOD  OF  PREDICTING  THE  INSTANTANEOUS  WIND  AT  THE  LAUNCHER  FROM 
SOME  REMOTE  POINT  IS  TO  FILTER  OUT  ALL  FREQUENCIES  WHICH  HAVE  NO  SIGNIFICANT 
CORRELATION  BETWEEN  THE  TWO  POINTS  WHILE  RETAINING  ALL  OF  THE  FREQUENCIES 
WITH  SIGNIFICANT  CORRELATION.  THE  COHERENCE  FUNCTION,  PREVIOUSLY  DISCUSSED 
IN  THIS  REPORT,  IS  A  MEASURE  OF  THE  CORRELATION  COEFFICIENT  FOR  ANY  PARTIC¬ 
ULAR  FREQUENCY.  WHEN  THE  COHERENCE  VALUES  ARE  KNOWN  FOR  THE  SEPARATION 
DISTANCE  BETWEEN  THE  SENSOR  AND  THE  LAUNCHER  THE  UNCORRELATED  HIGH  FREQUEN¬ 
CIES  CAN  THEN  BE  FILTERED  OUT  OF  THE  DATA  BY  USE  OF  A  MOVING  AVERAGE  OF 

PER  I OD  T  • 

C 

A  MOVING  AVERAGE  WILL  FILTER  OUT  THE  TURBULENCE  DUE  TO  FREQUENCIES 
HIGHER  THAN  1/(2^).  FOR  EXAMPLE,  FROM  FIGURE  5  AND  TABLE  IV  WE  SEE  THAT* 
UNDER  IDEAL  DAYTIME  CONDITIONS,  AN  AVERAGING  TIME  OF  ABOUT  26  SECONDS  IS 
NECESSARY  TO  OBTAIN  CORRELATIONS  OF  GREATER  THAN  0.7  FOR  A  HORIZONTAL 

(Cont.  on  Page  21) 
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Rocket  Wavelength 
Eddy  Wavelength 


ADMITTANCE  FUNCTION. 

Figure  6 
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Fig.  7 


SEPARATION  DISTANCE  OF  5^  FEET.  TABLES  I  AND  II  INDICATE  THAT  AN  AVERAGING 
TIME  OF  5  TO  10  MINUTES  IS  NECESSARY  TO  GET  A  RELIABLE  ESTIMATE  OF  THE  WIND 
FOR  A  HORIZONTAL  SEPARATION  DISTANCE  OF  1,000  FEET.  THE  COHERENCE  FUNCTION 
IS  ALSO  APPLICABLE  TO  ESTIMATE  THE  INSTANTANEOUS  WIND  IN  A  VERTICAL  PATH. 

Tables  I  and  II  indicate  that  an  averaging  time  of  approximately  5  minutes 
at  150  FEET  IS  NECESSARY  TO  PREDICT  THE  WIND  SPEED  AT  300  FEET. 


CONCLUS  IONS 


It  is  felt  that  low  level  WIND  PREDICTION  SCHEMES  FOR  LAUNCH  TIME 
ESTIMATES  OF  IMPACT  DISPERSION  PATTERNS  CAN  BE  IMPROVED  OVER  CURRENT  METHODS 
BY  USE  OF  THE  SPECTRAL  AND  CROSS-SPECTRAL  ANALYSIS  OF  THE  ATMOSPHERE. 

Simple  and  practical  methods  of  predicting  the  turbulence  field  for  field 

USE  HAVE  BEEN  DEVELOPED  FOR  PARTICULAR  LOCATIONS  AND  CAN  BE  ADAPTED  FOR 
ANY  LAUNCH  SITE  IN  QUESTION. 
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